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ABSTRACT 

We present sever al div erse applications of the spherical fast convolution method suggested by 
IWand elt & Gorski^ (|200ll), which is useful for studies of telescope optical properties and for con- 
struction of shaped filters for analysis of all-sky data. We study sidelobe pickup in three design 
concepts for the CALISTO infrared telescope. The beam convolution allows for direct comparison of 
the performance of each telescope design. At 100 microns, the best of these designs limits sidelobe 
contamination from Galactic dust emission to < 0.1 MJy/sr for most of the sky with \h\ > 25°. With 
the fast convolution method, we illustrate the impact of asymmetric primary beams on the recovery 
of the temperature power spectrum for the Planck microwave background mission. Finally, we use 
the fast convolution method to specify a class of orientable filters on the sphere, working through a 
pedagogical example. 

Subject headings: telescopes, methods: data analysis, infrared: general, cosmic microwave background 
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1. INTRODUCTION 

Observational astronomy and cosmology experiments 
require realistic simulation and analysis pipelines for de- 
tailed study of systematic effects. In particular, the sim- 
ulation of realistic telescope optical response is a compu- 
tationally expensive yet important task. Distortions of 
the beam (or point-spread-function) caused by imperfec- 
tions in the mirrors and other optical components of the 
telescope lead to undesirable systematic effects, includ- 
ing extraneous pickup by the beam far-sidelobes. These 
effects are important and need to be modeled. An exper- 
iment may observe a region far from the Galactic plane 
in order to avoid foreground contamination. However, 
scattered or diffracted Galactic emission may still spoil 
the measurement. At the same time, large numbers of 
past and planned large sky area surveys (2MASS, NVSS, 
SDSS, WMAP, Planck, LSST, etc.) emphasize the im- 
portance of analysis on the sphere. 

In order to fully simulate and analyze these effects, it is 
necessary to compute the detector response in every di- 
rection: convolving a 47r-steradian beam model with the 
sky for every possible pointing and orientation of the tele- 
scope. The computational cost of a brute force analysis 
in real space prohibits it for all but the poorest angular 
resolution. Alternatively, simplifying assumptions — like 
beam axisymmetry — cut the cost of harmonic space cal- 
culation, but lose the orientation information, a poten- 
tially significant failing^ 

Wandelt & Gorski (20011 ) devised a fast, efficient, and 
computationally feasible method for performing this type 
of all-sky convolution. This paper highlights this algo- 
rithm with a series of applications which are difficult or 
impossible to address by other means. In section [21 we 
describe the convolution method, and the preparation of 
the spherical harmonic representation of the beam and 
sky map. Although the applications in this work deal 
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with total intensit y only, this meth od has been extended 
to polarization by IChallinor et all (^000). 

This convolution method is used widely in two ways. 
First, as described above, it can simulate telescope opti- 
cal properties. It is particularly well suited to the study 
of large-scale far-sidelobe pickup. If the beam is not too 
small and can be described by a limited number of az- 
imuthal harmonics, this method is also suitable for sim- 
ulating the response of an asymmetric main beam. The 
simulation pipeline for the upcoming Planck Cosmic Mi- 
crowave Background satellite routinely uses it for this 
purpose (|Reinecke et al.|[2QQ6l ). In section [3l we offer two 
demonstrations of fast convolution probing optical prop- 
erties, describing the far-sidelobe pickup for the proposed 
infrared satellite mission CALISTO and the impact of 
main beam asymmetries on the temperature power spec- 
trum evaluation of Planck. 

Beyond telescope optics, the second main use for this 
convolution method is in the analysis of all-sky astro- 
nomical signals. For example, it has been used to 
l ook for Bianchi Vllh cos m ological mode l contributions 
(jJaffe et al.l I2005L l2006allbl : ICavon et all I2006D to the 
cosmic microwave background temperature map mea- 
sured by the Wilkinson Microwav e Anisotropy Probe 
{WMAP, e.g. iHins haw et al.l l2007f ). Also, in a num- 
ber of works, the method has been used to process all- 
sky signals with orientable filters, in preparation for 
examination of the CMB for no n- Gaussian signatures , 
cross-correlation, o r other studie s (iBarreiro et all 120061: 
Q-uz et al. 2006; Liu & Zhang "2006: "McEwen et al.l 
2006 , 2007a,b; Vielva et al. 2006a.b, 2007; Wiaux e tld] 
20061 . l2007al lbl). Often in these works, the convolution 



method is used to transport around the sphere a class 
of filters derived from projected wavelets. We discuss 
orientable filters on the sphere in section [3l where we 
emphasize that the convolution method is independent 
of the choice of (band-limited) filter; it does not depend 
on these wavelet-derived, or any other, particular filters. 
To make this point, we hide a non-trivial signal under a 
large amount of noise, and find it with an appropriate 
orientable, matched filter. 
Finally, discussion and conclusions follow in section [H 
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2. METHOD 
2.1. '^Total convolution^' 
The convolution problem evaluates 

T(n,c^) = J dn'[D{h,co)b]''{h') s{h'). (1) 

In this notation the vector n defines a position on the unit 
sphere defined by longitude (j) and colatitude 6^ and uj is 
the orientation at that position. The rotation operator 
D rotates the function 6(n) (the "beam") into position 
from a fiducial orientation at the North pole via Euler 
angles 0, cj), while the integral over the function s{h) 
(the "sig nal" ) covers the whole sph ere. For band-limited 
functions IWandelt fc Gorskil (|2QQli ^ found a class of fast 
algorithms for this type of problem, computing the con- 
volution over part or all of the sphere for some or all 
beam orientations. These employ the spherical harmonic 
representation of the beam and signal, defined by 



b{h) = ^ ^ bimYim{n), 



(2) 



^=0 m=-l 



for the beam (similarly for the signal) where L defines the 
band limit. For the special case considering the whole 
sphere at all orientations, dubbed "total convolution," 
the main result is 



T(n,C^)= ^rr 



^m\(j)-\-m29+m^(jO 



where 



mim2'm.3 



(3) 



(4) 



^=0 



Here A 



2 ~ ^mim2('^/2) reduced Wigner ro- 

tation matrices^ evaluated at it/ 2. For a uniform rec- 
tilinear grid in {(j)^9^uo)^ equation ([3]) may be computed 
rapidly using a fast Fourier transform. The main compu- 
tational cost is in computing T^im2m3 7 which scales as 
O(L^), compared to a computational cost of 0{L^) for 
a brute force approach. If the azimuthal content of the 
beam is limited, the scaling can be reduced by limiting 
ma, yielding a scaling of 0(I/^m3^max)- In the applica- 
tions in this work, the convolutions each require a few 
minutes to compute on modest computer hardware. By 
conventional means, these computations would each take 
several days. 

2.2. Preparation of the harmonic representation 

Because the convolution may be computed efficiently in 
a band-limited harmonic representation, we need to con- 
struct a harmonic representation of the beam (and sky) 
which captures the relevant aspects of the function on 
the sphere, but does not introduce unphysical features, 
such as spikes at the sampling positions of the beam or 
ringing due to a sharp band-limit imposed in harmonic 
space. 

Often the beam and sky maps are presented as equally 
spaced samples on isolatitudinal rings. Compared to an 



^ See lEdmondsl (| 19601, c hapter 4) for the definition of the d- 
functions and IRisboT (|1996l ) for recursion relations for rapid com- 
putation. 



arbitrary pixelization, this arrangement makes harmonic 
analysis more efficient: it allows the azimuthal portion 
of the transform to be computed by an FFT and re- 
duces the numbe r of Leg;endre function evaluations (see 
e.g. HEALPix,^ Gorski et al."^2005). Additionally the 
samples may be distributed azimuthally on meridians. 
Because physical optics calculations of the telescope ra- 
diation pattern are resource-intensive, a practical exam- 
ination of the sidelobes may sparsely sample the sphere. 
Sidelobes typically subtend large solid angles, so this is 
not usually an issue, but it does have some implications 
for the evaluation of the beam harmonics. 

The resuperposition of the harmonics (possibly on a 
different pixelization) will attempt to reproduce the orig- 
inal sampling of the beam as a collection of (5-functions. 
Due to the finite band limit, these ^-functions are some- 
what smeared. However, if the sampling density of the 
output convolution anywhere exceeds that of the input 
grid, these sampling spikes will be visible. This common 
problem arises when the sampling is sparse, or when the 
areal density of samples varies widely, such as on any grid 
on which the number of samples per latitude ring is con- 
stant, like the geographic grid. Smoothing the harmonics 
can help address these issues with the beam representa- 
tion, and yield a result closer to the physical beam. 

A useful smoothing for grids is an azimuth-direction- 
only smoothing: 



d<p'f{4>-4>')biO,<P) 



(5) 



where the smoothing kernel has J^^ d(j) f {(/)) = 1. If we 
expand the kernel as /(</>) = /m exp(— zm(/)), the 
spherical harmonic coefficients of the smoothed function 
are. 



bin 



d{cosO)d^b{e,^)YC^{e,^) 



= 27Tfmblm' (6) 

This result is straightforward, and may be obtained by 
substitution of the definition of 6, expanding / and b in 
harmonic series, integrating over the azimuthal angles, 
using the orthogonality relation of associated Legendre 
polynomials, and summing over the resulting Kronecker 
^-functions. 

One particularly useful example smoothing kernel is 
the unit-integral boxcar. 



m 7^ 0. 



(7) 



(8) 



{2M~\ -(/>o < < (/>o < ^; 

0, otherwise, 
for which the harmonic coefficients are 

{1, m = 0; 
sin(m0o) 

Another useful smoothing function is a symmetric 
Gaussian. After smoothing with a symmetric Gaussian, 
the spherical harmonic coefficients are 

km flkm (9) 

where 

/,=exp(-Za + l)cTLooth/2) (10) 
Note that neither of these filters modify the beam 
monopole. Thus the integral under the beam is unaf- 
fected by the smoothing. 

^ http://healpix.jpl.nasa.gov/ 
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TABLE 1 

Beam parameters 





^max (°) 


AO (°) 




<^smooth ( ) 


on- axis 


1 


5.0 X 10-4 


0.0833 


5.0 X 10-4 




25 


8.0 X 10-4 


0.0833 


1.0 X 10-3 




180 


5.0 X 10-4 


360.0 


2.0 


off- axis 


3 


5.0 X 10-4 


0.0833 


5.0 X 10-4 




40 


1.0 X 10-3 


0.0833 


1.0 X 10-3 




180 


1.4 X 10-3 


5.0 


2.0 



Note. — Pixelization of the beam for on- and off- axis 
CALISTO design cases from the GRASP computation. Each 
is divided into three non-overlapping latitude zones, marked 
by the maximum colatitude (^max) in each zone. The pixel 
centers in each zone are separated by and Acf). In prepar- 
ing the convolutions, the beams in each region were smoothed 
by an azimuthal boxcar filter of width A(f) and a symmetric 
Gaussian of width cTsmooth- 

3. APPLICATIONS 
3.1. Far-sidelobe pickup: CALISTO 

For our first application of the above convolution 
method, we examine telescope far-sidelobe pickup. In 
particular, we studied design concepts for a space-based 
infrared telescope. The Cryogenic Aperture Large In- 
frared Space Telescope Observatory, or CALISTO^ is 
tuned for the wavelength range 30-300 fim. The whole 
optical system, including the ~ 5 m diameter primary 
mirror, is cooled to 4 K, and yields a main beam reso- 
lution of ~ 1.2''-12'', depending on wavelength. Based 
on the low emission of the cold optics, a very low sys- 
tem noi se equivalent power can be achieved, ^ 10~^^ W 
(jGoldsmith et al.ll2006l : iGoldsmith et al.ll2007[ ). 

The purpose of this study was to calculate and com- 
pare the sidelobe pickup for three optical designs: (1) an 
on-axis design, where the secondary is held in front of the 
primary by three struts, (2) an off-axis Gregorian design, 
and (3) an off-axis Gregorian design with a 4 K cold stop 
placed between the detectors and the secondary. The 
sidelobe pickup is a convolution of the sidelobes with 
the sky, and is a straightforward application of equation 
([!]). The all-sky beams were calculated previously us- 
ingGRASP9 in three latitudinal zones (| Goldsmith et al.l 
|2QQ7[ ) (see table [T]), with decreasing angular resolution 
away from the main beam, which is centered at ^ = 0°. 
The on-axis beam contained ~ 139 million pixels, and 
the off-axis beams contained ~ 193 million pixels. 

In figure [TJ we show the beams from each of the three 
designs. Each design shows spillover from the primary 
mirror (where the spillover is opposite the axis of the 
main beam) and the secondary (where it is closer to the 
main beam). In the on-axis design, the support struts 
cause a 3-fold symmetric reflection pattern. The off-axis 
designs have lower sidelobes, particularly in the case in- 
cluding a cold stop. 

To suppress the main beam, in each case we multiplied 
the inner 20 arcminutes by an azimuthally symmetric 
Gaussian weight function, smoothly rising to unity at 
6>o = 20' with a standard deviation of ag = 2': 




Fig. 1. — The beams for three optical designs compared in the 
CALISTO sidelobe study. To the left, the beams are displayed 
in the Gnomonic projection, with the main beam at the center. 
To the right, the beams are shown in the Mollweide projection, 
with the main beam at the top. Color scale indicates log^oC^) and 
saturates in the sidelobes of the on-axis design. (Top) The on- 
axis design, with notable reflections. (Middle) The off-axis design 
with no cold stop, where the spillover from the secondary (closer 
to the main beam) and primary are not aligned along the optical 
axis. (Bottom) The off-axis design with a cold stop is similar, with 
reduced sidelobes. 

The harmonic components were computed separately 
in the three latitude zones. We wrote a custom code 
to compute the spherical harmonic transforms, using a 
harmonic tran sform library de veloped within the Planck 
collaboration ([Reineckd l2006l ) and MPI^ for multipro- 
cessor computation. We experimented with tapering the 
transition between the latitudinal zones for spherical har- 
monic transformation, but found this was not required to 
capture the large scale features of the beam. 

In the spherical harmonic basis, the three zones were 
then separately smoothed with a pixel-scale azimuthal 
boxcar and a symmetric Gaussian, and then summed. 
(See table [U for smoothing parameters.) The smoothing 
is necessary to suppress spikes in the resuperposition of 
the beam, caused simply by the input pixelization (from 
the GRASP computation). Figure [2] show the result of 
this smoothing. Without it, the spikes around the cen- 
ters of the original pixels propagate into the convolution, 
leading to unphysical structure in the result. 

To perform the convolution, we use the 

IWandelt Gor ski' (20011) ^^totalconvol yer" software 
in its implementation optimized by iReinecke et al.l 
(|2006l . section 6.5). We limit the harmonic content to 
L < 512, although we limit the azimuthal beam index 
ms < 256. On a 4-processor, 2.4 GHz Opteron machine, 
the convolution completes in 3 minutes. We repixelize 
the convolution into a set of N^ide = 512 (7' pixel) 
HEALPix maps, one for each azimuthal {co) rotation of 
the beam. 

For visualization, we take 2-dimensional slices from the 

^ http:/ /www-unix. mcs.anl.gov/mpi/ 



Fig. 2. — Resuperposition of harmonic components for the off-axis, cold stop case. Color indicates \ogiQ{b). Values less than zero are 
shown in gray, and are an unphysical consequence of the finite band limit, but are small enough to have little consequence in the convolution. 
(Left) Without smoothing; note the shape of the low-resolution sidelobe, like a 5-function comb. (Right) The smoothed representation. 




Fig. 3. — The sky and the sidelobe pickup for prospective CALISTO designs. Color indicates logiQ(MJy/sr). F or the conv olutions, slices 
through the convolved sky at fixed orientations are depicted as spheres in the Mollweide projection. (Top) Schl egel et al.l (1998) galaxy 
model at 100 microns. (Middle-left) A slice through the convolution for the on-axis design. (Middle-right) A slice for the off-axis design 
with no cold stop. (Bottom-left) The off-axis design with a cold stop, with the spacecraft oriented to minimize sidelobe contamination in 
the northern galactic hemisphere. (Bottom right) The same, but oriented to minimization contamination in the southern hemisphere. 
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3-dimensional convolution in {(j)^6^uj). At fixed cj, we 
plot the convolution as a map on the surface of the sphere 
defined by {4>^0). That is, each map pixel displays the 
entire contribution to the instrument response from the 
beam sidelobes, when the beam center is directed at the 
pixel and rotated about that axis according to uo. For 
a (5-function signal on the = plane, maps of these 
slices look like images of the beam, pinwheeling as the 
orientation angle is incremented. This can be misleading, 
causing the viewer to think the map is a composite image 
of beams instead of a slice of the convolution. This dis- 
tinction becomes clearer near the pole, where the coordi- 
nate singularity of this representation distorts the beam 
image. For example, the convolution of any beam with 
a ^-function at the north pole will appear azimuthally 
symmetric in these slices with fixed uo. 

In Figure [31 we sh ow sample convolution s with the 100 
micron dust map of ISchlegel et al.l (|l998l ). The on-axis 
case shows the most sidelobe pickup. More importantly, 
rotating the instrument about the optical axis has lit- 
tle impact on the structure of the sidelobe pickup. The 
off- axis cases, however, give much better control over the 
placement of the sidelobes. Because the sidelobes are not 
symmetric about the optical axis, there is the opportu- 
nity to orient them away from the galaxy. For most of the 
sky with galactic latitude \h\> 25°, the cold- stop design 
allows an orientation which limits the foreground contri- 
bution to < 0.1 MJy sr~^, or roughly a few hundredths 
of a mJy per main beam. 

3.2. Scanning with asymmetric main beams: Planck 

For a scanning experiment with asymmetric beams, the 
instrument response for every observation can be chal- 
lenging to compute. The resulting maps depend both 
on the beam shape, and the details of the scanning path 
and orientation, which can serve to symmetrize the ef- 
fective beam. Because the fast convolution method si- 
multaneously computes the beam response in all pos- 
sible configurations, it provides an attractive approach. 
Here the challenge is that for a small main beam, the 
band limit must be high to fully capture the beam struc- 
ture. However, since the main beam often has limited 
azimuthal content (see lFosalba et all (|2002l ) for the case 
of an elliptical Gaussian), even with a high band limit, 
the technique of section [2] may still be tractable. This ap- 
proach is routine ly u sed in the Planck simulation pipeline 
(jReinecke et al.ll2006f ) to incorporate the effects of asym- 
metric beams. 

Here we quantify the impact of beam asymmetries on 
one of the main observables for Planck., the temperature 
power spectrum. We use simulated data from the 30 
GHz channel of the Low Frequency Instrument (LFI), 
wi th 4 detectors at 30 G Hz (the simulations are detailed 
in lAshdown et a/]l2007l ). These beams have a full-width 
half maximum of ~ 33^ and have substantial asymmetry, 
since the horns are ne ar the edge of the focal plane (the 
Planck "Bluebook," iThe Pla nck Cohaborationl 120061) . 
The beams are modeled as elliptical Gaussians; as an 
indication of realistic beam asymmetry, we t ake an axis 
ratio similar to the K band from WMAP (|Page et al.l 
[2003). We compare to axisymmetric Gaussian beams 
with the same effective area. 

Noiseless time-ordered data were generated with the 
convolutions of the beams and the CMB sky, interpolat- 




FlG. 4. — The impact on the power spectrum due to asymmetric 
beams. (Bottom) the power spectrum of the CMB, with no beam 
apphed, after simulated scanning with the four detectors of the 
Planck LFI 30 GHz channel with symmetric Gaussian beams with 
FWHMs of ~ 32^, and scanned with asymmetric elliptical Gaus- 
sian beams with the same effective area. (Top) the ratio of the 
asymmetric to symmetric beam transfer functions, corresponding 
to the square root of the power spectrum ratio. 

ing the response of each detector along the 1-dimensional 
scan path through the 3-dimensional convolution. The 
realistic scan pattern had the satellite spin axis pointing 
along a cycloidal path around the anti-sun direction for 
one year. The convolution band limit was L = 3000 
and the maximum azimuthal index of the beam was 
^3, max = 14. This procedure is handled by software in 
the Planck simulation pipeline. Here, the time-ordered 
data are simply binned to make a pixelized map of the 
sky (HEALPix N^ide = 512), then the power spectrum is 
evaluated. 

Figure [H compares power spectra. The asymmet- 
ric beams are compared to symmetric Gaussian beams 
scanned the same way. At very large scales, the beam 
and its shape are immaterial. At intermediate scales, up 
to / ~ 400, the asymmetric beam suppresses power more 
than the symmetric one. By contrast, at small scales, the 
symmetric beam suppresses more power, since the short 
axis of the asymmetric beam is able to sample smaller 
scales. 

3.3. Orientable filters 

For our third application, we examine fast convolu- 
tion as a method to search for specified patterns in data 
on the sphere. In sever al works which process the sky 
with steerable wavelets (Ba rreiro et al.|[2Q06l: ICruz et al.l 
20061: iLui &: Zhang 2006; McE wen et al.l l2006l \2007M : 



Vielva et al. 2006a.b, 2007; Wia ux et al.ll200 



the~'Wandelt & Gorski (2001) convolution mechanism 
underlies the orientable filtering, making these meth- 
ods computationally efficient and numericall y tractable. 
The sp irit of the formalism here is similar to iJaffe et al] 
(|2006q ). 
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Suppose our data set is a scalar function d{n) on the 
sphere, the sum of the signal and Gaussian noise with 
noise covariance C(n,n'). For the case that we know 
the shape of a signal template, so(n), and we want to 
estimate the amplitude {A) of this template, our model 
for the signal is A[Dso]{h^), where D = D(n, cj) is the 
rotation operator. For computational ease, we briefly 
switch to a vector notation, viewing our data as d{h) = 
d. The likelihood is given by 

- 2 log£(A|d, n, cj) = (d - A[L)so])^C-i 

{d-A[Dso]) (12) 

plus an immaterial constant. We can write the maximum 
likelihood estimate of A: 



A{d^ n, co) 



[Z)so]tC-i[Dso] 



(13) 



Viewed as a filter acting on d to produce A, this is known 
as a matched filter. 

Switching back from our vector notation, the numera- 
tor is easily written as a convolution (equation [T]) : 



/ 



= J dh' [D{h,co)sonn){C-'d){h'). (14) 

If the noise is statistically homogeneous and isotropic, 
the integration over the inverse covariance is most easily 
done in spherical harmonic space. 



(c-M) 



Im 



dim 

Ci ■ 



(15) 



Equation ([Hj) may now be computed using the total con- 
volution method. 
The denominator of equation [13] is 



dh'dn'' [Dsonn')C-\h\h'')[Dso]in'') (16) 



Again if the noise is isotropic, we may drop the com- 
mon D rotations of the signal template and evaluate the 
denominator in harmonic space, where it is equal to 



;=o m=-i ' ;=o 



(17) 



In figures [5] and [H we illustrate an example of the ori- 
entable matched filter. We have hidden a beam with a 
nontrivial shape under white noise. Only at the largest 
harmonic scales does the signal exceed the noise. For 
this reason, a simple large Gaussian filter can discover 
the general location of the signal, but orientable match- 
ing filter finds it with much better precision, and can 
identify the orientation to high significance. 

4. DISCUSSION 

We have explored three applic ations of the sp herical 
fast convolution formalism of Wandelt fc Gorskil ([2001). 



First, we examined the sidelobe pickup from CALISTO, 
quantifying the level of contamination in three optical 
designs from galactic duct emission at 100 microns. We 
found that the off-axis design with a cold stop permitted 
only a low level of contamination away from the galactic 
plane. In future we wish to address the issue of side- 
lobe pickup from zodiacal light. We note that the side- 
lobe techniques described here are not at all specific to 
C A LIS TO, and may be applied to any experiment. In 
particular we are interested in further applications of this 
technique to the upcoming Planck mission when detailed 
sidelobe beam maps become available. 

Second, we examined the impact of asymmetric beams 
on scanning experiments, showing the difference in the 
smoothing of the CMB power spectrum as measured by 
Planck from scanning with symmetric or asymmetric ide- 
alized beams. Characterizing the beam asymmetries, 
and the impact on the maps, power spectrum, and as- 
sociated errors, looks to be one of the most challenging 
data analysis problems for the Planck mission. The type 
of work displayed here is only a first step in approach- 
ing this problem. Correcting this beam asymmetry effect 
remains an open research issue. 

Third, we looked at orientable matched filters, showing 
that the use of orientation information can greatly aid 
the recovery of specific templates from noisy data. We 
believe these techniques can find use in characterizing 
directional asymmetries and to look for systematics for 
all- sky surveys. 

Finally, we believe that the fast convolution formalism 
is — and will continue to be — a useful tool for many data 
analysis applications in spherical geometry. At this time, 
an optimized software implementation is being readied 
for public release. 
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Fig. 6. — Signal and noise power spectra for the case in figure \E\ (top-left), and slices through the convolution fixing all angles but 
one. These slices pass through the maximum of the convolution, and show the recovery of the position and orientation of the signal at 
((/),^,cj) = (270°, 120°, 210°). Red is for signal alone; green includes noise. 
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